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Abstract 

By demanding a compact spectrum for the right-handed neutrinos and an 
, approximate quark-lepton symmetry inspired from SO{W) gauge unification 

(assuming a Dirac neutrino mass matrix close to the up quark mass matrix), 
we construct a fine tuning scenario for baryogenesis via leptogenesis. We 
find two solutions with a normal hierarchy, with the lightest neutrino mass 
mi different from zero, providing an absolute scale for the spectrum. In the 
approximations of the model, there are three independent CP phases : 5l 
(that we take of the order of the quark Kobayashi-Maskawa phase) and the 
^ two light neutrino Major ana phases a and /3. A main conclusion is that, 

^ although this general scheme is rather flexible, in some regions of parameter 

c5 space we find that the necessary baryogenesis with its sign is given in terms of 

the 5l phase alone. The light Majorana phases can also be computed and turn 
out to be close of 7r/2 or small. Moreover, 50(10) breaks down to the Pati- 
Salam group 5'C/(4) x SU{2) x SU{2) at the expected natural intermediate 
scale of about 10^'^ — 10^^ GeV. A prediction is done for the effective mass in 
(/3/3)o!/ decay, the mass and the sum of all light neutrino masses. 
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1 Introduction and qualitative remarks 



The discovery of oscillations, advocated so many years ago by Pontecorvo [T], in 
solar and atmospheric neutrinos is one of the most important experimental discov- 
eries of the last century, the most relevant after the proposal of the Standard Model 
and its precision tests. The discovery of neutrino oscillations is also a milestone in 
the search of New Physics (NP). 

Up to now four quantities related to the Pontecorvo, Maki, Nagakawa and Sakata 
(PMNS) matrix [2] [3] have been experimentally measured : 

Ami - 8 X 10"^ (1) 

tan^ 9s ~ 0.4 (2) 

Ami - 2-5 X 10"^ eV^ (3) 

tan^ Oa ^ 1 (4) 

where the subindices s and a mean respectively solar and atmospheric neutrinos. 
An upper bound has been been found for the component of Uei along the heaviest 
niass eigenstate 

sin^ ^13 < 0.05 (5) 

and the limits 

m^^ < 2.2 eV (6) 
\ <mee>\ < 0.4 eV (7) 
E < 1 (8) 

i 

from the high energy spectrum of the electrons in nuclear beta decay, from the upper 
limit on the rate in neutrinoless double beta decay (for Majorana neutrinos) and 
from astrophysics. 

Interestingly, a more restrictive bound combining all cosmological data has been 
obtained recently by G. Fogli et al. |1] : 

Y,m,^<0.2eV (9) 
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to which we will refer in Section 8, comparing it to our results. 

But for the moment, in this qualitative introduction, we will rely on the generally 



The most natural framework to account for the order of magnitude of neutrino 
masses is the seesaw model [3] , where the 6x6 neutrino mass matrix has the form 



where the 3x3 Dirac neutrino mass matrix m^i has elements of the order of the 
masses of charged fermions and is the Majorana mass matrix of the right-handed 
neutrinos, which are singlets of the Standard Model gauge group, with elements of 
the order of the scale of breaking of the lepton quantum number. 

The information on oscillations gives us only four of the nine parameters of 
the light neutrino mass matrix. Within the simplifying assumption of neglecting 
6'i3 and consequently the neutrino Dirac CP violating phase, we will be able to 
strongly constrain the value of its smallest eigenvalue, and fix the values of the two 
higher ones, as well as the two Majorana phases, simply by demanding that these 
parameters have a soft dependence on the values of the matrix elements of M^. 
We will obtain these results despite the fact that we expect a rather hierarchical 
spectrum for the eigenvalues of itld, as it happens for the other fermions and as is 
natural in a 5*0(10) framework. The mathematical principle is quite simple : it is 
that the inverse of a function with a critical dependence on a variable is a very slowly 
varying function : the product of the derivatives is of 0(1). The demand of having 
matrix elements and eigenvalues of Mr of the same order, given a mixing matrix of 
leptons similar to the one for quarks, will fix mi and the Majorana phases of light 
neutrinos. As a result of this requirement, we shall get a compact spectrum for the 
Nr masses, which will make the leptogenesis scenario for baryogenesis natural, as 
well as predictions for the electron neutrino mass bounded from tritium (5 decay 
and the matrix element | < m^e > \ appearing in neutrinoless double beta decay. 
By compact spectrum for the heavy right-handed neutrinos we simply mean to have 
eigenvalues of the same order of magnitude. 

From the seesaw formula 



accepted loser bound 




(10) 
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rriL 



(11) 



one gets 



det Mr 



(det mo) 



det iriL 

From eqn. dsj) one obtains the upper limit 



(12) 



det mi I < — eV^ 



(13) 



while in principle there is no lower hmit for the Lh.s. of the inequahty (13). Notice 



that we write the absolute value in the l.h.s. of (13) because neutrino masses, being 



Majorana masses, can differ in sign for neutrinos with opposite CP. 
Moreover, from eqns. ([T]) and ^ we get : 



Am! 



R2I 



\mi\ 



8 X 10"^ eV 



Ami 



\ms\'-cos% \m,\'-sinX ^ 



2.5 X 10-'^ eV^ 



(14) 
(15) 



where the unfamiliar formula (15) for Am^, proposed in [6], is demonstrated in 



the Appendix. This formula is an improvement over the usual ones found in the 



literature, Am^ 



|m2p or Am^ 



R3I 



I mi p. Of course, in the limit 



\mi\ all these formulas coincide. However, we must underline that the results 



of this paper are not really sensitive to adopting formula (15) or the usual ones 



From the preceding formula one gets a lower limit for the ratio 



|m2 I 
ms 



> 0.18 



(16) 



A temptative lower bound for | detmj;,| may be found in the 50(10) framework 
by taking, as in [7j, 

Idetm^l = 4 X 10-2 Gel/3 (^^^^^ 
and for | det Mr\ the upper limit 



det Mr\ < 2.7 X 10^^ GeV^ 
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which comes by assuming that the three right-handed neutrinos take a mass at the 
scale of B — L spontaneous symmetry breaking in the 5*0 (10) model, with breaking 
to the SU{A) X SU{2) x SU{2) Pati-Salam group [8] at the intermediate scale 3 x 10^^ 

Emni. 



We then get, from the seesaw formula (12) : 

Idetm^l > 6 X 10"^^ e\/=^ (19) 
Assuming a normal hierarchy for the light neutrinos : 



mal ~ J Ami - 8-9 x 10~^ eV (20) 



mal ~ J Ami - 5-0 x 10"^ eV (21) 



eqn. (19) will then imply the following lower bound for \mi\ : 

|mi| > 1.3 X 10"^ e\/ (22) 

i.e., a non-vanishing value for the lightest neutrino mass mi, an absolute scale for 
the light neutrino spectrum. 

As we will see below, a rather sharp prediction for mi and relevant predictions 
for the l.h.s. of eqns. (|6])-([8]) will be achieved by our demand of a compact Mr 
spectrum and successful leptogenesis. 

The measurements of the cosmic microwave background anisotropies jTT] and 
the abundance of light nuclei produced in primordial nucleosynthesis [12] give a 
consistent value for the baryon asymmetry : 

Yb = — ~ ~ 9 X 10"^^ 23 

s 7.04 ^ ' 

This baryonic asymmetry may arise from the leptogenesis scenario [13], with 

a leptonic asymmetry produced at a high scale, which gives rise by the B — L 

conserving sphaleron processes [H] at the electroweak scale to a baryon asymmetry 

below that scale. 

Within the leptogenesis scenario, the baryon asymmetry, baryon to entropy frac- 
tion, is given by 



5 



Yb ^ -\yl (24) 

that should be compared with the experimental value given by (23). 

Concerning Grand Unification, the SU{5) minimal model is disfavored, since it 
generates a small baryon asymmetry at the high scale, washed out at the electroweak 
scale, since in that model B — L is conserved. Thus, 5*0(10) with its 5 — L generator 
spontaneously broken, that we will adopt in its non-Supersymmetric version, should 
be preferred to SU (5) to realize the leptogenesis scenario. 

The paper is organized as follows. In Section 2 we give the relevant formulas 
for the inverse seesaw, mass matrices and mixings. In Section 3 we formulate our 
SO{10) Ansatz. In Section 4 we give the formulas needed for CP violation and 
the baryon asymmetry. Section 5 is devoted to a simple mathematical procedure 
to obtain a quasi-degenerate right-handed neutrino spectrum (that presents a level 
crossing) and a realistic light neutrino spectrum. We underline an illuminating limit 
of considering, for the matrix diagonalizing rriD, a. pure Cabibbo matrix that then 
we extend to a general matrix of the CKM form, introducing therefore CP violation. 
We find two possible solutions. In Section 6 we expose a simple procedure to slowly 
lift the degeneracy of the heavy right handed neutrinos, and give the corresponding 
evolution of Am^ and Am^. In Section 7 we exhibit the results for CP violation 
and baryon asymmetry, in the one-fiavor approximation, and in Section 8 we give 
the predictions for rriu^ and the effective neutrino mass in {f3f3)ou- In Section 9 we 
relax a reality assumption used in Sections 6 and 7. In Section 10 we comment on 
the compact heavy neutrino spectrum and on the level crossing region. Finally in 
Sections 11 and 12 we underline open problems within the present approach and we 
conclude. 



2 Inverse seesaw, mass matrices and mixings 



From (11), we can deduce the inverse seesaw formula. 



Mr = -m 



D 



m, 



mo 



(25) 
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and diagonalizing the neutrino Dirac mass matrix m^i by 



(26) 



one gets the formula [H] 



Mr = -m'j, ml^ mD = - m^"' l^^* m^^ \/^+ m^''' 
where the last equality follows from the definition [B] of the matrix 



(27) 



= V^* ml^ 

The neutrino mass matrix rriL is diagonalized by the PMNS matrix U 



(28) 



where 



and U writes : 



rriL = U* mf'"'^ 



diag(mi,m2,m3) 



(29) 



(30) 



U 



C12C13 S12C13 Si3e 

"■S12C23 — Ci2S23Si3e*'' C12C23 — Si2S23Si3e*'' S23C13 

i5 „ r. 0^0 ^*<5 



\ S12S23 - Ci2C23Si3e'" -C12S23 - Si2C23Si3e'" C23C13 / 

. diag(l,e^",e^^) (31) 

where 5 is the Dirac phase and a and /3 are the Majorana phases. For the latter we 
adopt the convention of Davidson et al. [13J. 

Let us make a remark on the counting of phases. One has, in all generality, 6 
independent phases in the Type I seesaw scheme, as established in [15] [I6] and as 
it is exposed in the review [13] (last reference. Section 2.1). In the model that we 
develop below, the number of independent phases will be reduced according to the 
hypotheses adopted. 

Taking into account the data on solar and atmosperic neutrinos and the fact 
that Si3 is bounded to be small, we will take 
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Sl3 ^ 

and approximate, from now on, the matrix U as follows : 

U - - % % 4^ . diag(l, e*", e'^) 



(32) 





Ss 





Ss 


Cs 


1 


V2 






Ss 




1 



(33) 



and CP violation originating in the Dirac phase 5 drops out. 

To simplify the expressions in what follows, we change the notation for the 



diagonal matrix in (29)-(31) 



diag(mi,e *°m2,e 



diag(mi,m2,m3) 



(34) 



where, from now on, m2 and are assumed to be complex parameters. 



Eqn. (29) now writes, in the approximation (33), and with the notation conven 



tion of the r.h.s. of (34) : 



rriL 







V2 



Cs 

V2 



\ V2 



1 
V2 

V2 V2 / 



diag(mi, 777,2, ms) 



Cs - 

Ss 





V2 V2 



V2 
1 
V2 



Cs 

72 ) 



(35) 



Let us strongly underline again that in (35), and in what follows, the parameters 



m2,ms are assumed to be complex, containing, according to (34), the Majorana 



phases defined by (31 ). 

These phases will be computed at different stages. Their calculation will depend 
on some hypotheses to be made explicit below, and on the successive parametriza- 



tions assumed for the matrix rriD (26) 



From the previous hypotheses we obtain the following complex symmetric matrix 



?77, 



Cs I -^s 

nil m2 



V 



CaSs ( A 

{---) 



) K 



1 

rrt2 
1 



V2 



CsSs I 1 

I Cg I 
mi r?i2 ma 

1 _L _£l L 

2 \mi m2 ms 



1 

mi 



V2 

_1 (i _L ^ _ 
2 \mi m2 

2 \mi m2 



1 

m2 

]_ 

ma 

1 

ma 



and the matrix Al ( 28 ) is also complex symmetric 



(36) 
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A 



Lt 



A' 



(37) 



In a previous work [6j, to comply with the lower bound for the mass of the lightest 
right-handed neutrino claimed by |17], we did set upper limits on the coefficients of 
contributions proportional to the products of the Dirac matrix eigenvalues (m^jg)^ 
and m£)2^Di in the Mr matrix, related to uil by the inverse seesaw formula. 



From formula (27) we see that to get a quasi-degenerate heavy Majorana neutrino 
spectrum we need that the terms proportional to (m^g)^ and m^igm^g have to be 
small, that means for a hierarchical Dirac mass spectrum that the matrix elements of 
(28) A33 and = have to be small. From (28), we get the following expression 



for these matrix elements of interest : 



2 \ mi 1712 rns 



\mi 




a/2 Vmi 



m2 



21 



1 
2 



mi 1112 



V^2* + 



2 V 

1 ^ "2 



1 

m2 

■^s _|_ 

nil ^2 



v.. 



23 



St ct 1 

+ — + — 



yL* 



2 \mi 777.2 



V, 



23 



(38) 



^33 



32 



31 



5_ + 

777i 7772 



^/2 \mi 7772 



a/2 \mi 1712 



yL* 
'^31 



^^31 



a/2 Vmi 7772 
2 \ 777i 7772 ^3 



^30* + 



C e S a 



^^32 



31 



1 
2 



7771 ^^2 



32 



a/2 Vr77i 

1 ^ "2 



1 

7772 
777 i 7772 



1/, 



33 



1 

2 V 7771 



— + — 
7772 "^3 



33 



33 



(39) 



3 Our 50(10) Ansatz 



The seesaw model is realized in the framework of 5*0(10) unified gauge theories 
|18j . where 5 — L is a generator, which has to be spontaneously broken. Long time 
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before the firm evidence for neutrino oscillations this phenomenon had been claimed 
|19j as the most promising experimental signal for 5*0(10) unification. 

A systematic study of the spontaneous symmetry breaking in 5*0(10) unified the- 
ories has lead to propose [9] the model with SU (4) x SU (2) x SU (2) [8] intermediate 
gauge group, broken at the scale of order 3 x 10^^ GeV [lOJ. 

A general analysis has been done in [7] on the possibility to construct a realistic 
leptogenesis scenario within the seesaw model with neutrino Dirac masses in a hier- 
archical ratio, as it is the case for u-type quarks. The most promising case has been 
found with M3 ~ 10^^ GeV and nearby values for the masses of the two lightest 
right-handed neutrinos. 

Although in the present paper we follow the general idea [7] of leptogenesis 
generated by quasi-degenerate right-handed neutrinos, we look for a more compact 
spectrum for Nji, with the heaviest right-handed neutrino at the intermediate scale, 
of the order 10^^ GeV. 

In 5*0(10) the hypothesis that the electroweak Higgs transforms as a combination 
of 10 representations implies at the unification scale the equalities among mass 
matrices 



rrie = rud (40) 
niD = niu (41) 

For b and r masses relation ( 40 ) at the intermediate scale is in reasonable agree- 
ment with experiment but, as Georgi and Jarlskog pO] have shown in the SU{5) 
case, one needs also higher dimensional representations. The generalization of this 
argument to 50(10) was given by Harvey et al. [21]. For an overview on fermion 
masses and mixings in gauge theories, see the review article [22] . 

Within 50(10), with the electroweak Higgs boson belonging to the 10 and/or 
126 representations, and no component along the 120 representation, the mass ma- 
trices are symmetric. As a consequence, the unitary matrices V^ and V^ that 



diagonalize Dirac neutrino matrix ( 26 ) are related : 



yR ^ yL* 



(42) 



and the matrix M/j (27) becomes 
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Mr = -mn mr^ 



yL+ ^diag j^L ^diag yL* 



(43) 



Let us now go back to the question of the phase counting, quoted for the seesaw 
scheme in Section 3, in the particular case of 5*0(10) with symmetric Dirac neutrino 



matrix. Since V has only one phase, and m^, through the mixing matrix U (31) 



has three phases, we have reduced the number of independent phases, from 6 in the 



general case to 4 independent phases. In the approximation (32) S13 ~ that we 
have adopted, this means that we have 3 independent phases, namely a phase from 
V^, that we will call 5li and the two Majorana phases a and /3 from (33). 



Below, in Sections 5 and 6, we will impose two other conditions that reduce 
further the number of independent phases, from 3 to a single one. 
For the diagonalized Dirac neutrino matrix 



/ 



diaq 



V 




















(44) 






mo, J 

we will adopt the numerical values proposed in [7] , inspired from the up-quark mass 
matrix : 



10-3 GeV 



tud^ = 0.4 GeV 



mo-. 



100 GeV 



The matrix m^"^ m'^"'^ appearing in (27) has the form 

/ 

diaq \ L diaq 

rrij-, A rrij^ = 



43 



rriD^mD^ A^^ mo^rnD^ A^ ^ 
y mDirriDs A^^ mo^mD^ ^23 



rriD-^mD^ A 



L 

12 



22 ^D2^Ds ^23 



^D-i ^33 



(45) 



(46) 



that clearly shows that in order to have a compact N^i spectrum from (27) one needs 



small values for the matrix elements and ^33. 

For V^ we will assume a form qualitatively similar to the Cabibbo-Kobayashi- 
Maskawa (CKM) quark matrix, that reads, in the standard convention (except for 
the phase 5^, we take the same notation as for the light neutrino mixing matrix 



(31), but in what follows there is no ambiguity) 
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Cl2Cl3 
~'5l2C23 ~ Ci2S23Si3e 



C12C23 ~ ■5l2'S23'Sl3e 



\ S12S23 - Ci2C23Si3e 

where Sl is the CP-violating phase. 



-C12S23 - Sl2C23Sl3e' 



■523C13 
C23C13 J 



(47) 



Formula (47) is correct at least for exact quark- lepton symmetry, with Higgs 



in the 10 representation : if the mass matrices (40) are diagonal and real one has 



= VcKM- For phenomeno logical purposes we assume this form in what follows. 



We define, as usual, in terms of Wolfenstein parameters : 



S12 = A 



S23 = AX'^ 



(48) 



Of course, in our problem the parameters X,A,p,ri do not necessarily have the same 
precise values as in the quark sector : we are interested only in an order of magnitude 
estimate. 

Let us say some words concerning the diagonalization of the right-handed neu- 



trino matrix. Since in our 5*0(10) Ansatz Mr (43) is complex and symmetric, we 
can diagonalize it by using a single unitary matrix : 



Mr = Wr M^^"'^ W'r (49) 

The matrix Wr is such that all eigenvalues are real and positive. The effect of 
phases will appear in the matrix Wr. These phases will of course have consequences 
for baryogenesis and for neutrinoless double beta decay. 

As we will see below, our demand of suppressed values for A^^ and ^23 generates 
a compact form for the Mr spectrum, which helps in getting in a natural way the 
desired lepton asymmetry. 



4 Leptogenesis and baryon asymmetry 

In this Section we recall the basic formulas concerning the CP violating asym- 
metry ei and the corresponding baryogenesis asymmetry Yb^ . We work in the basis 
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in which the mass matrices of charged leptons and of right-handed neutrinos are 
diagonal, i.e. from (25) and (49) : 



M' 



diag 
R 



m* m-' run 



R 



(50) 



Therefore, in the computation of the CP-violating asymmetry ei we define 



rriD = mo 



(51) 



such that 



By convention we label the masses of the heavy neutrinos Nr. (i = 1, 2, 3) : 



< Ml < Ma < Ms 



(53) 



In terms of niD, the CP asymmetry writes, for the lightest heavy neutrino Nr-^^ : 



1 



E/ 



M2\ Im \{m^mD)lk 



(54) 



where v = 174 GeV is the scale of electroweak symmetry breaking, and the function 
f{x) is given by [23] : 



f{x) = y/x 



1 — X 

that in the limit a; >> 1 becomes : 



1 , X , /I + a; 

+ 1 - 1 + X lo! 



X 



(55) 



2x/x 



fix) 



and the effective neutrino mass, that controls the amount of washout, writes 



(56) 



mi 



Ml 



(57) 



The cases that we encounter in our calculations below satisfy the strong washout 
condition 



mi » 3 X 10"^ eV 



{51 



13 



and the corresponding baryon asymmetry writes, in the one-flavor approximation, 
that we will adopt in the following [21] : 

. J 0.3 ^1 ^'""'^ I (59) 

2 fi-* V "^1 / 
where g^, ~ 107 in the Standard Model, in the non-Supersymmetric case. 



5 Quasi-degenerate heavy right-handed neutrinos 
and realistic light neutrino spectrum 



In order to get a compact A^^^; spectrum, a sufficient condition is to impose that 
the matrix elements A33 and are suppressed, because we are dealing with the 



matrix (46) and the niD eigenvalues (45). As a first exercise, we thus consider the 



solutions of the equations, linear and homogeneous in the inverse of the neutrino 
masses ^ (^ = 1,2,3), 

A23(mi, ma, m-i) = A^^{mi, mg, mg) = (60) 

We are aware that this is a very drastic assumption, but will help to guide our 
research of a compact right-handed neutrino spectrum, and also to look for its 
consequences on the light neutrino masses and the amount of baryogenesis that 
one can get. We must emphasize that in this Section, and in the following ones, 
we are dealing with a fine tuning scheme. We cannot content ourselves with just 
order-of-magnitude estimates, but we need precise numerical calculations. 



Notice a new important point in the phase counting of eqn. (43) with the hy- 



pothesis (32). Under the two reality conditions (60) that we now impose, the 3 
phases (see Section 3) are now reduced to a single phase, either 6l or one of the two 
Majorana phases a or /3. 

Since we do not have experimental information on the Majorana phases, we 
will, from now on, compute a and /3, and later the CP asymmetry ei and baryon 
asymmetry Yb^ in terms of 6l- Of course, in principle one could also compute the 
pair {6l, a) in terms of (3 or P) in terms of a. But in the present 5*0(10) 
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approach the natural thing to do is to take as input 5l, since we can take it to be of 
the order of Kobayashi-Maskawa (KM) phase 5kmi on which we have information. 

5.1 in the limit of a pure Cabibbo matrix 



For our purpose, it is a good illustration to study the consequences of this hypothesis 
considering it within the very simplified approximation of a 2 x 2 Cabibbo matrix 



^ ^ 



-Sl2 Ci2 

1 



(61) 



Since from (61) 5l drops out, we are left with only two phases, namely the 



Majorana phases a and (3. Imposing the two reality conditions (60), these phases 



will be fixed, as we see below. 



From (60) and (61), we find that the light and heavy neutrino spectra turn out 
to be reasonable. From (38) (39), the matrix elements of we are interested in are 



A 



23 



A 



A 



33 



32 - 

1 
2 



1 
2 



mi m2 



C e S e 



Cl2 



\mi 



1 



1 

m2 



Sl2 



mi 1712 ?7T,3 



(62) 



If we impose the very strong assumption (60) we have the two equations 

-,2 



mi m2 



mi m2 



tan^ 



12 







S's Cl 1 

— + — + 

mi m2 ms 



(63) 



and solving for m2 and m^ in terms of mi, one gets : 



m2 



\/2 — tan^s tan 6 



12 



a/2 taxi^Os + tan6's tan 6 



mi 



12 



ma 



V2 



iaxiOs tan 6* 



12 



mi 



(64) 

idJiOs tan ^12 

From the data ([l])-(|4]) and formulas (14)(15) one gets, from (64), a number of solu- 
tions for mi and 6*12. 
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However, if one looks for solutions with 612 in the neighborhood of the Cabibbo 
angle 6c, one finds the following solutions, according to the sign of tan^^, since only 
its square ^ is measured : 

(1) For tan ~ — -^0.4 one gets. 



mi 



0.0030 eV 



7712 



tan 012 = 0.140 
= -3.1522 mi 



(65) 

ms = -16.9273 mi (66) 

We have taken several digits for the mj values to get a compact spectrum for 
the Nr^ and, as we will see below, for later to obtain the non-degeneracy of the 
two higher states. The reason is that we are dealing with a fine-tuning problem. 
Of course, one could take the first digits and say the degree of approximation at 
each stage, but we believe that our way of presenting the results, although harder 
to read, corresponds better to the reality of the calculation. 



In consistency with (34) we assume the convention mi > and one gets the 
following spectrum 



mi = 0.0030 eV 



m2 = -0.0094 eV 



= -0.0507 eV (67) 



where two heavier neutrinos have opposite CP from the lighter one. This means 
that the Majorana phases are 



TT 



a 



f3 



TT 



(68) 



We find, from (26), (42) and the value tanO^ (65), for the symmetric Dirac mass 
matrix : 



rriD 



0.0087 



-0.0549 



-0.0549 0.3923 
100 



GeV 



(69) 



and 



Ml = 5.5504x10^ Gel/ M2 = 1.42991 x 10^° GeV" M3 = 1.42992 x 10^° Gel/ 

(70) 



(2) For tan^. 



'OA one gets 
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tan5i2 = 0.243 (71) 
|mi| = 0.0062 el/ ma = -1.752 mi mg = 8.191 mi (72) 

Assuming again the convention mi > one gets the following hierarchical spectrum 



mi = 0.0062 eV 



m2 = -0.0109 eV 



mg = 0.0509 eV (73) 



where two neutrinos (the lightest and the heaviest) have opposite CP from the third 
one. This means that the Majorana phases are 



a 



TT 

2 



13 = 



(74) 



We obtain for this solution, from the value tan ^12 (71), the Dirac neutrino matrix 



/ 



0.0233 



-0.0916 



V 



GeV 



-0.0916 0.3777 
100 

and the quasi-degenerate right-handed heavy neutrino spectrum : 



(75) 



Ml = 6.72168 X 10^ GeV M2 = 8.30366 x 10^ GeV M3 = 8.30409 x 10^ GeV 

(76) 

The results for these two solutions seem encouraging because we get in both cases a 
value of the angle ^12 that is rather close to the Cabibbo angle 6c- It seems highly 
non-trivial and amazing that such a simplified form of the V^ matrix could give 
already these results consistent with quark-lepton symmetry. 



5.2 with approximate CKM form 



We now switch on the other V parameters and consider the full matrix (47). To 



perform the calculations we adopt for mi and tan 6*12 the values obtained in the 
pure Cabibbo limit, namely mi = 0.0030 eV, tan ^12 = 0.140 for solution (1) and 
mi = 0.0062 eV, tan ^12 = 0.243 for solution (2). 

For the numerical calculations we thus proceed in the following way : 
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i) From either solution (1) or solution (2) of the preceding Subsection, we fix the 
parameters mi and tan 612 : 



tan 9^ 



^0.4 



mi = 0.0030 eV 



tan^i2 = 0.140 (77) 



(2) tan ^, = +7(14 mi = 0.0062 eV" tan^i2 = 0.243 (78) 

ii) For the rest of the matrix elements we deduce the Wolfenstein parameter 



A in (48) from (77) and (78) and take, just a guess, the parameters A, p and rj from 



the quark sector CKM matrix, i.e. for example 



A = O.i 



p = 0.13 



T] = 0.35 



(79) 



and, using these parameters, we fix S23, S13 following (48) 



;i) tan^ia = 0.140, S23 = 0.0154, S13 = 0.0008, 5^ = 1.2152 (80) 



gives 



( 0.990 
-0.139 
\^ (18.6 -7.4i) X 10" 



0.139 (2.8 
0.990 
-0.015 



- 7.5i) X 10-^ ^ 
0.015 

1- ) 



^1) 



while from 



(2) tan 012 = 0.243, S23 = 0.0446, S13 = 0.0039, (5l = 1.2152 (82) 



one obtains 



0.972 
-0.236 

\^ (9.20 - 3.58i) X 10-3 -0.044 



0.236 (1.37 - 3.68i) X 10-3 \ 

0.971 0.045 
3 _n n/i/i \ 



13) 



iii) Then, we solve equations (60) for m2 and m3 and compare with experiment 
for Ami and AmJ,. 



Notice that this numerical procedure is less rigid that the one adopted in the 
simpler case of a pure Cabibbo matrix of the preceding Subsection, where Am^ 
and Am^ were fixed to the experimental central values ((l])(|3| and we did solve for 
tan^^i2, mi, m2 and m3. We prefer to change here our numerical approach due to 
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the extreme fine tuning of the problem. It seems to us sensible enough if we get 
results for Am^ and Am^ that are roughly consistent with experiment. 



We find the following results. 



(1) For tan 6*^ ^ -VOA, mi = 0.0030 eV and tan 6*12 = 0.140, one gets : 



ma = -0.0095 e°-°°^^^ eV 
that correspond to the Majorana phases 



-0.0495 e 



0.0075i 



eV 



a = ^- 0.0018 



13 = - -0.003^ 



^5) 



Let us notice an important point. We obtain the CP violating part of the Majo- 



rana phases for the light neutrinos (i.e. their departure relatively to | in (85)) from 



the 6l phase, that we take close to the KM phase 6km- This can seem paradoxical, 
because 6^ concerns the Dirac neutrino mass. However, because of 6l, the matrix 
is complex. This implies that, setting mi real as we have done above, the so- 



lutions from eqns. (60) for m2 and ms (with the notation (34)) must be complex. 
The departure of these Majorana phases relatively to the ones obtained in the real 



Cabibbo limit (68) (74) turn out to be numerically small. 



We obtain, from (14) (15) and (84) : 



Ami = 8-1 X 10"^ eV^ Ami = 2-4 x 10"^ eV^ (86) 

The agreement with the data is good. 

Let us now give the complex symmetric Dirac neutrino mass and the right-handed 
heavy neutrino spectrum for this solution. We get 



0.0090 + 0.0003i -0.0576 - O.OOlli 0.1849 + 0.0739i 



mD= -0.0576 - O.OOlli 0.4155 - 0.0003i -1.5206 + 0.0103i GeV{87) 
0.1849 + 0.0739i -1.5206 + 0.0103i 99.9764 
and the quasi-degenerate right-handed heavy neutrino spectrum : 

Ml = 5.53144x10^ GeV M2 = 1.43230x10^° GeV M3 = 1.43232x10^° GeV 
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(2) For tan 61^ ^ +V0^, mi = 0.0062 eV and tan6'i2 = 0.243 one gets : 



m2 = -0.0106 e"° °^^^ eV = 0.0455 e° °°^^' eV 

that correspond to the Majorana phases 



(89) 



a = -| + 0.0080 



We obtain, from (14) (15) and (89) 



13 = -0.0039 



(90) 



Ami = 7.4 X 10"^ eV^ 



Ami 



2.0 X 10"^ eV^ 



(91) 



The agreement with the data is not as good as for solution (1). We could change 
the initial conditions for mi and tan and get a better agreement. However, it is 
not our intention to make a fit but to get a qualitative agreement with the data. 
We get the Dirac matrix for this solution 



mo 



( 



\ 



0.0304 + 0.0066i -0.1319 - 0.0148i 0.9152 + 0.3575i 
-0.1319 -0.0148i 0.5676 - 0.0076i -4.3450 + 0.0869i 
0.9152 + 0.3575i -4.3450 + 0.0869i 99.8003 



GeV (92) 



and the quasi-degenerate right-handed heavy neutrino spectrum : 



Ml = 6.84678 x 10^ GeV M2 = 8.84878 x 10^ GeV M3 = 8.84909 x 10^ GeV 

(93) 



The signs and phases of the results (84) and (89) will have quantitative consequences 



for the effective neutrino mass in neutrinoless double beta decay, as we will see below. 



5.3 CP violation and baryon asymmetry 



The results of the preceding Subsection show that imposing the very drastic condi- 



tions (60) one gets quasi-degenerate right-handed neutrino spectra. 



To have a feeling on how to proceed, let us make an exercise in the case (1), where 



the quasi- degeneracy (88) is less pronounced. Let us compute the Dirac matrix (51) 



in the basis in which the heavy right-handed neutrino mass matrix is diagonal, ei, mi 
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and finally Yg^ . We will assume that the lightest neutrino decays out-of-equilibrium 
and that one can apply the one-flavor approximation. 

We find the following result 



mo 



-0.055Z -0.052 + 0.132Z -0.131 - 0.052i 

-0.001 + 0.391i -0.006 - l.OSli 1.080 - 0.006i 
y -0.002 + 0.347i -0.064 + 70.702i -70.702 - 0.064i y 



GeV (94) 



Of course, unlike expression (87), this matrix is no longer symmetric. Using it in 



eqns. (54), (57) and (59), we obtain 
ei ~ -3.805 X 10"^° mi '. 



0.050 eV 



3.05 X 10 



-15 



(95) 



where we have used the exact formula (55) for the function f[x) since the three 



heavy neutrinos are rather close in mass, although Nr^ is lighter, and formula (59) 



is applied because, according to the value of mi (95), we are in the strong washout 



regime (58). 



The obtained baryon asymmetry Yb^ = 3x 10^^^ is much too small, by about four 
to five orders of magnitude, although of the right sign. The reason is the smallnes 
of ei, that follows from the quasi-degeneracy of M2 and M3 and the opposite CP 
asymmetry contribution from both heavy neutrinos. Indeed one finds, for the two 



terms in (54) : 



/ 



) 87rf2 (mi^m^) 11 



-/ 



Ml) Sttv^ (m^mi) 11 



-3.634 X 10 



-6 



(96) 



that shows a strong cancellation giving a very small CP violation ei. 

Although for the moment we do not get good phenomenological results, we should 
however emphasize an interesting limit of the present scheme, namely : 

Sl implies ei Ybj 



a = j3 — (solution (1)) 



a ^ , /3 -> (solution (2)) (97) 



Notice that we have an intuitive argument to understand the quasi-degeneracy 



between M2 and M3 (88) and the smallness of the CP violation (95). Indeed, in the 
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limit = = (60), and neglecting terms of order m|,^ and m£,^m£,^, since we 
take to be close to a diagonal matrix, the Mr matrix has only M^^g, Mr,^^ and 
Mi^22 sizeable matrix elements with Mr^^ ^ ^Rsii ^-nd the matrix Mr is close to 
real. Therefore one has M2 — M3 and ei ~ 0. 

We can try to modify our very simplified scheme by lifting the degeneracy of 



Nr^ and Nr^. We will thus relax somewhat the strong condition (60), but keeping 
the general physical idea of a compact heavy Nr spectrum. We will allow for non- 
vanishing values for the matrix elements 1^23 1 ^^"^ I ^33 1 5 keeping them "small", i.e. 
values much smaller than each of their individual contributions that, due to the 
smallness of the light neutrino masses, are naturally of the order ~ 10^^ GeV~^, as 



can be seen in eqns. (38) (39) 



As we will further examine, one can thus obtain a rather compact Nr spec- 
trum, and also reasonable values for the baryon asymmetry consistent with the data 
without spoiling the good properties of the light neutrino spectrum. 

6 Lifting the quasi-degeneracy of heavy neutrinos 



We will proceed now, in terms of some parameters, to a continuous and slow 
lifting of the quasi-degeneracy of the heavy right-handed neutrino masses obtained 



in the previous Section within the strong hypothesis (60). 

In this Section we do the calculation considering non-vanishing values for the 



r.h.s. of the eqns. (60). Moreover, since we have seen that even the drastic as- 
sumption of taking ^33 = A23 = gives reasonable neutrino spectra (84) or (89), 



we will allow to vary IA23I and IA33I within a very wide range, keeping "small" 
values (<< 10^^ GeV~^), and observe how the heavy neutrino and the light neu- 
trino spectra evolve, as well as the consequences for the CP violation asymmetry 
ei, the effective neutrino mass rhi and baryon asymmetry Yr^- We will perform the 
calculations for both solutions (1) (tan c:^ —\/QA) and (2) (tan6's ~ +\/()A). 

are linear homogeneous 



Notice that the strong conditions (60) A 



L 

33 



^23 



equations in 



1,2,3). We now allow for no n- vanishing inhomogeneous terms 
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A^^{mi,m2,m3) = C^3 (98) 
Af3(mi,m2,m3) = ^3^3 (99) 

To lift the very close degeneracy between M2 and M3 in the case examined before, 
A33 = ^23 = 0, we just need to have non- vanishing, in general complex parameters 



C23 and Cj"3 in the r.h.s. of (98) and (99). However, to have an overall compact 



heavy neutrino spectrum we need inhomogeneous terms that should be small in 
modulus relatively to each individual term in A33 and A23 (38)-(39). This means 



that we will take non-vanishing values for C^g and C^^^ with the condition 

1^2^31. 1^2^31 « 10^^ GeV-^ (100) 
In principle one should scan the general two complex numbers C23 and C^^ and 
see how the heavy neutrino spectrum evolves, as well as the light neutrino masses, 
the light neutrino Majorana phases a and /3, the CP asymmetry ei and final baryon 
asymmetry Yg^. 



As pointed out at the beginning of Section 5, the Majorana phases (85) and 
(90), that we found for C23 = C^a = have their origin in the approximation 



adopted for the matrix V^, that we take close to the CKM matrix. In order to 
preserve this interesting feature, we will assume that the non-vanishing values of 



the inhomogeneous terms C^3 and C^;^ are real and satisfy (100). Later, in Section 
9 we will relax this reality assumption and will see that we have a wide domain of 
values for complex Cjjj and C^^ that can give reasonable results. 

The first important observation to be made is that the degeneracy between Nfj,^ 
and Nfi.^, that we have found solving eqns. (60) is lifted considerably if jC^sl is non 
vanishing in some region with |C^3| << 10^^ GeV~^, and we have realized that the 
mass difference M3 — M2 is rather insensitive to the precise value of |Cj"3|, provided 
that its value is not "too large". Importantly, the amount of CP violation, and 
therefore of baryon asymmetry depends also on the value adopted for |C^3|. 

To reduce the number of parameters, we assume that C^3 and Cjjj are real and 
of equal modulus |C^| = |C;^3|. As we vary and C^^g we find essentially the same 
heavy neutrino spectrum and roughly the same values for Ami ^^"^ ^'^a? indepen- 
dently of their relative sign. We find that what is dependent on this relative sign is 



23 



the amount of CP violation and baryon asymmetry. If C^g = C^^, (independently of 
its sign), the baryon asymmetry can be at most of 0(10^^^), but if C^g and C^g are 
of opposite sign, one can get a correct amount of baryon asymmetry. 

In conclusion, after some trial and error guesses, we respectively adopt real num- 
bers for and C^g for both solutions (1) (tan6's < 0) and (2) (tan^^^ > 0) 



(1) - = Cg^g > (101) 

(2) - C4 = Cg^g < (102) 

with ICj'gl << 10^^ GeV . As it will become clear below, the adopted sign for each 
of the solutions corresponds to the experimental sign > in some region for the 
parameters C^, C;fg. We assume that Yb^ — Yb, an hypothesis that will be justified 
in Section 10. 

We will now show how the heavy neutrino and the light neutrino spectra evolve 



under the conditions (101)(102). In the next Section we will show how ei, rhi and 
behave. 

Of course, our ansatz for C^g, C^g is just a guess. We do not intend to make a fit 
to the overall data, light neutrino spectrum and baryon asymmetry. We just want 
to see if the description of these data is possible within this scheme of a compact 
heavy neutrino spectrum and approximate quark-lepton symmetry. Other equations 



of the type (98) (99), with complex r.h.s. values for the parameters Cj^g, Cj"g give also 
acceptable results, as we will see in Section 9. 



6.1 Heavy and light neutrino spectra for case (1) tan^^ < 

To perform the calculations we adopt again the values of the pure Cabibbo limit, 
namely mi = 0.0030 eV, tan 6*12 = 0.140, although one could slightly change these 
values to get a better fit. 

In Fig. 1-a the right-handed heavy neutrino spectrum is plotted. In Figures 1-b 
and 1-c we show respectively the solar and atmospheric quantities Am^, Am^. 



24 



14 
12 
10 
8 




8 



10 



Fig. 1-a. Log-log plot of the right-handed heavy neutrino spectrum (masses in 
GeV units) as a function of — C^^g = > in units of GeV"^, for fixed mi = 
0.0030 eV and tan = 0.140. Within the range -C^^ = Cj^g = 10^ - 10^ GeV-'^ 
there is a level crossing. The angular points come from the fact that the curves are 
obtained from interpolation of a finite number of points. The same applies to the 
other figures. 
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Fig. 1-b. Ami units as a function of = C33 > in units of GeV ^, 

for fixed mi = 0.0030 eV and tan 6*12 = 0.140, in a log scale for C^^- 
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Fig. 1-c. Ami units as a function of — = C^^ > in units of GeV ^, 



for fixed mi = 0.0030 eV and tan6'i2 = 0.140, in a log scale for C. 
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Let us comment on these figures. The angular points that appear in the figures 
are an artifact of the representation of the curves, obtained from an interpolation of 
a finite number of points. Notice that for each point we must perform the singular 
value decomposition of the matrix Mr in order to compute the quantities necessary 
to obtain the baryon asymmetry. 



The first striking point is that, as we have learned from eqns. (98) (99), the A''^ 
spectrum (Fig. 1-a) is very compact for C^^ not "too large", C^fg < 10^ GeV"^. As 
explained in the Introduction, there is an expected correlation between the stability 
of the light neutrino spectrum and the compact heavy right-handed neutrino one. 
The fine-tuning for the close heavy neutrino masses ensures the stability of the light 
neutrino ones. For — = C'gg > 10^ GeV~^ the right-handed neutrino spectrum 
evolves into a hierarchical spectrum. The values obtained for Am^ (Fig. 1-b) and 
Am^ (Fig. 1-c) are very stable and consistent with experiment for a wide range of 
values of — = C^^, of about eight orders of magnitude. 

We observe two other important things in Fig. 1-a : the degeneracy between Nr^ 
and Nji^ is lifted, and there is a level crossing around — ~ C^^ = 3 x 10^ GeV~^. 

An important point to be also underlined is that one of the levels [Nji^ before the 
crossing) remains practically constant in the whole studied range, while the mass of 
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Nji^ decreases. After the crossing we will call Npii this right-handed neutrino, being 



the lightest, according to convention (53). 

As we can see from Figs. 1-b and 1-c, the values obtained for Ami ^^"^ 
in good agreement with the data for a very wide range of the parameters. 



For the Majorana phases a and /3, shown in Figs. 1-d and 1-e we find rather 
constant values (in a logarithmic scale) that are very close but a little smaller than 
|, as shown in the figures below. 
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Fig. 1-d. The Majorana phase a as a function of — C^3 = Cjjj > in units of 
GeV-\ in a log scale for Cg, for fixed mi = 0.0030 eV and tan ^12 = 0.140. The 
X-axis is centered at 7r/2. 
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Fig. 1-e. The Majorana phase ^ as a function of — C^3 = C^3 > in units of 

GeV~^, in a log scale for C^g, for fixed mi = 0.0030 eV and tan6'i2 = 0.140. The 
X-axis is centered at tt/2. 

6.2 Heavy and light neutrino spectra for case (2) tan 6*5 > 

To perform the calculations we adopt again the values obtained in the pure Cabibbo 
limit, namely mi = 0.0062 eV, t&nOu = 0.243. 

In Fig. 2-a the right-handed heavy neutrino spectrum is plotted. In Figures 2-b 
and 2-c we show respectively the solar and atmospheric quantities Am^, Am^. Figs. 
2-d and 2-e display the result for the Majorana phases a and (3. 
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Fig. 2-a. Log-log plot of the right-handed heavy neutrino spectrum (masses in 
GeV units) as a function of = — > in units of GeV~^, for fixed mi = 
0.0062 eV and tan ^12 = 0.243. Within the range -C^j = 10^ - 10^ GeV'^ there is 
a level crossing. 
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Fig. 2-b. Am^ in eV"^ units as a function of C^g = —C^^ > in units of GeV~^, 
in a log scale for —C^^, for fixed mi = 0.0062 eV and tan 6*12 = 0.243. 
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Fig. 2-c. Am^ in eV"^ units as a function of C^g = —C^^ > in units of GeV ^, 
in a log scale for —C^^, for fixed mi = 0.0062 eV and tan^^i2 = 0.243. 

As shown in the figures below, in a logarithmic scale, we find for the Majorana 
phase a a rather constant value that is very close but a little larger than — | and 
for /3 a small negative almost constant value. 
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Fig. 2-d. The Majorana phase a as a function of C. 
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— Cjj, > in units of 



GeV-\ in a log scale for -Cg^g, for fixed mi = 0.0062 eV and tan ^12 = 0.243. The 
X-axis is centered at — 7r/2. 
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Fig. 2-e. The Majorana phase /3 as a function of C^g = 
GeV-\ in a log scale for -Cg, for fixed mi = 0.0062 eV and tan ^12 = 0.243 



The first striking point in this case is that, as imposed from eqns. (98) (99), 
the Nfi spectrum (Fig. 2-a) has the same features as for the precedent solution, 
although it is very compact for —6*^3 < 10^ GeV^^, much more than in case (1). 
Within the range C23 = —^33 = 10^ — 10^ GeV~^ there is also a level crossing, 
on which we will comment below. For = — Cj^ > 10^ GeV~^ the right-handed 
neutrino spectrum evolves also into a hierarchical spectrum. Secondly, Am^ (Fig. 
2-b) and Am^ (Fig- 2-c) are very stable for a wide range of values of C^3 = —C^;^, 
of about seven order of magnitude. However, the agreement with experiment is not 
as good as for solution (1), although it is acceptable within a 3a range. Of course, 
we could somewhat change the initial conditions mi = 0.0062 eV, tan0i2 = 0.243 
and become closer to the data. This could be done, but we will not do it because 
our purpose is only a qualitative one within our (fine-tuning) scheme. 
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7 CP violation and baryon asymmetry in the re- 
gion approaching the level crossing 



Let us turn now to the quantities that are important for Baryogenesis via Lep- 
togenesis. Labelhng the hghtest heavy right-handed neutrino Nfi-^, our calculations 
show that the quantities — ei, rhi and Yb^ have a strong discontinuity across the 
level crossing region. We call always Nn^ the lightest heavy neutrino, even after the 



crossing, according to the convention (53) 



We will justify and characterize this term of level crossing, and discuss its impli- 
cations before this region in Section 10. 

To simplify the presentation of the results, we will restrict ourselves to the region 
before the crossing, where M2 and Mi become relatively close, i.e., to the following 
regions, slightly different in both cases : 

(1) 10^ GeV-^ < -C4 = eg < 10^-^ GeV-^ 

10^ GeV < M2 - Ml < 8.3 x 10^ GeV (103) 



(2) 10^ GeV-^ < eg = -eg < 10^-^ GeV-^ 

10^ GeV < M2 - Ml < 2. X 10^ GeV (104) 

To avoid the delicate situation related to the quasi-degeneracy of two heavy 
neutrinos, extensively studied by A. Pilaftsis et al. p5], we need the condition (see 
also [7]) 

Ti « M2 - Ml (105) 

where Fi is the width of the lightest heavy right-handed neutrino, that has an upper 
bound qualitatively given by [7j : 

r. < l^M. (106) 

Before the level crossing region one gets, from the parameters quoted above {rrit ~ 
rriD, ^ 100 GeV, v = 174 GeV and Mi ~ 5. x 10^ GeV) : 
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Ti < 3. X 10^ GeV (107) 
Therefore, before the level crossing region, taking into account the inequalities 



(103) (104), we see that in both cases (1) and (2) the condition (105) is satisfied. We 
are far away from resonant leptogenesis and we do not have to face complications 
related to quasi-degenerate heavy neutrinos, for which M2 — Mi < Ti. 



Let us now show the quantities — ei, mi and Yb^ for both solutions within the 



interesting ranges (103)(104). 



7.1 Case (1) tan^^ < 



Fig. 1-d displays the CP violation parameter — ei. Fig. 1-e the washout parameter 
mi, and in Fig. 1-f the baryon asymmetry Yb^ in the one-flavor approximation. 
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Fig. 1-e. rhi in eV units as a function of — = C33 in units of GeV~^, in a log 
scale for Cg, for fixed mi = 0.0030 eV and tan ^12 = 0.140. 




We observe that ei is negative and becomes large enough in absolute magnitude 
to give a large positive with rather stable values of the parameter rhi that imply 
strong wash-out in the whole region. Notice that — ei as well as Ybi grow as the 
mass difference M3 — M2 slowly grows and the mass difference M2 — Mi becomes 
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smaller, a phenomenon already underlined by Akhmedov et al. [7]. 

The main conclusion that we can draw from Fig. 1-f is that there is no problem 
to get a baryon asymmetry of the right order of magnitude {YB)exp — 9 x 10^^^. 



We must underline that if we did took the opposite sign in (101), we would 
have obtained the opposite sign for ei and therefore also for Yg^. Therefore, our 
scheme does not predict the sign of Yb since it depends on the chosen sign of the 
inhomogeneous terms. 



7.2 Case (2) tan^^ > 

Fig. 2-d displays the CP violation parameter — ei. Fig. 2-e the washout parameter 
rhi, and in Fig. 2-f the baryon asymmetry Yb^ in the one-flavor approximation. 




35 



1.2 


1 


1 


/ 


0.8 




. 6 


J 


0.4 




0.2 





4.25 4.5 4.75 5 5.25 5.5 



Fig. 2-e. mi in eV units as a function of Cj^j = — C^3 in units of GeF ^, in a log 
scale for -C^.^, for fixed mi = 0.0062 eV and tan ^12 = 0.243. 
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Fig. 2-f. Log-log plot of Yb^ as a function of = —C^z in units of GeV ^, for 
fixed mi = 0.0062 eV and tan Q^i = 0.243. 



We observe that e\ is negative and becomes large in absolute magnitude to give a 
rather large positive with values of the parameter mi within the strong wash-out 
regime. 

We must again point out that if we did took the opposite signe of the r.h.s. of 
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(102), we would have obtained the opposite sign for ei and therefore also for Yg^. 

The main conclusion that we can draw from Fig. 2-f is that in this case we are 
somewhat short of having a baryon asymmetry Yb^ of the right order of magnitude 
O^B)exp — 9 X 10^^^. However, as pointed out above, one could modify the initial 
conditions (the values mi = 0.0062 and tan^^i2 = 0.243) and get results in better 
agreement with the data. But it is not our purpose to make a detailed fit for Am^, 
Am^ and Yg^, we want just to give a qualitative trend. 

Let us emphasize again that in the present scheme developped in Sections 6 and 
7, due to the reality conditions on and C^^ we have again the interesting limit 



(197|) : 

Sl ^ implies ei -> Yb^ 
a = /3^^ (solution (1)) a -> , /3 (solution (2)) (108) 



8 Results for m,j^, (P/3)ou and sum of neutrino masses 

We give here the predictions for the electron neutrino mass, on which one has 
limits from tritium P decay : 

TTT-ue = cos^ 6's \mi\ + sin^ 6's \m2\ (109) 

and for the effective mass relevant for neutrinoless double beta decay | < rriee > \ 



that writes, within the approximation (33), 



< rriee > I = I cos^ mi + sin^ 6's m2 1 (HO) 



The neutrino masses and their phases for both solutions (84) (89) are very close 



to those obtained in the region that give an acceptable value for Yb- Taking thus 



the values for both solutions, using the notation (34) : 

(1) mi = 0.0030 eV, ma = -0.0095 e°-°°^^^ eV, = -0.0495 e°°°^^^ eV (111) 

(2) mi = 0.0062 eV, ma = -0.0106 e-° °^^' eV, mg = 0.0455 e°-°°^^^ eV (112) 
we obtain, respectively : 

(1) m^^ ~ 4.9 X 10"^ eV | < mee > | ~ 5.7 x 10"^ eV (113) 

37 



(2) m^^ ~ 7.5 X 10^^ eV | < m^e > | ~ 1.4 x 10"^ eV (114) 

For both solutions, due to the relative signs among the (i = 1, 2, 3) (i.e., due 
to the Majorana phases), there is a strong cancellation between the two terms in 



(110), a phenomenon already exhibited in |6j in another context. The cancellation 
is stronger for solution (1). 

For the sum of the absolute magnitude of all neutrino masses, we obtain : 

(1) 5] |m^J = 0.0620 e\/ (115) 

i 

(2) ^ |m^J = 0.0623 el/ (116) 

i 

One gets very close results for both solutions, that comply with the cosmological 
bounds (|8|) and (|9|) [4|. 

Let us make a last qualitative remark comparing the different possible future 
experiments on neutrino masses and stress the importance of cosmological limits. 

If one takes mi ^ one finds, from the data, |m2| — wAm^ ~ 9 x 10~^ eV 



and Imal ^ y Am^ + cos^^sAm^ ~ 5 x 10~^ eV, which correspond to a value for 
the l.h.s. in eqns. ([8])([9]) of the order 6 x 10~^ eV, very near the value that we 
have found and only a factor 3.3 below the bound ^ [4J- So, according to the 
present scenario, the most promising search for effects of neutrino masses, apart 
from oscillation experiments, is the analysis of cosmological data, while for beta 
decay and neutrinoless double beta decay one should need an improvement of more 
than two orders of magnitude. 

9 Relaxing the additional reality constraints of 
the model 



We now relax the conditions of the particular model that we have studied quan- 
titatively, namely given by eqs. (glQ with real and Cg^g satisfying (p^([l02|, 
and allow complex numbers for these parameters, keeping however "small" values 



for the moduli, as stated in (100). 
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Notice the important point that now we have two new sources of CP violation 
besides a single independent phase 5l (or a or and we recover the situation in 
which there are three independent phases : 5l, a and /3. However, as we will see 
below, in the range of interest for Yg^, the magnitude of the Majorana phases is not 
very different than in the real case studied in detail in Sections 5, 6 and 7. However, 
their new contributions have important implications for the baryon asymmetry. 

Just to have a feeling of what can happen, we take an extreme case and adopt 



relations (101) (102) with the condition (100), but taking now C^g and Cj'g purely 



imaginary. Interestingly, the results are phenomenologically good and show that 
our general scheme of a compact Nr spectrum is flexible enough. 

We do not give the corresponding curves of Sections 6 and 7, and give values for 
representative points with acceptable phenomenological results. 

For case (1), i.e. tan ~ -^04 with mi = 0.0030 eV and tan 6*12 = 0.140, 
taking 

- C^, = C^3 = ^ 10-^ (117) 

we find the following results : 

Ami = 8-1 X 10"^ eV^ Ami = 2-4 x 10"^ eV'^ 

a=- -0.0018 /3 = - - 0.0038 

2 ^2 

Ml = 5.531 X 10^ Gel/ M2 = 1.383 x 10^° Ge\/ M3 = 1.483 x 10^° Ge\/ 
mi = 0.050 eV ei = -2.755 x 10"^ Yb, = 2.211 x 10"^° (118) 



while for case (2), i.e. tan9s ~ +VOA with mi = 0.0062 eV and tan 6^12 = 0.243, 
taking 

- C^3 = C^3 = 10-' (119) 

we find : 



Ami = 7.4 X 10"^ eV^ 



Ami = 2-0 X 10~3 eV^ 



a = -- + 0.0079 



(3 = -0.0039 
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Ml = 6.847 X 10^ GeV M2 = 8.844 x 10^ GeV M3 = 8.954 x 10^ GeV 



rhi = 0.170 eV 



-3.622 X 10^ 



Yb^ = 7.035 X 10 



-11 



(120) 



These results are phenomenologically reasonable, and we find a whole region in 
their neighborhood that gives also good results. 

Notice that in the numbers of case (2) we are not far away from saturating the 



bound (107), and therefore we are approaching the regime of resonant leptogenesis. 

Let us emphasize again that in the case examined here we have two different 
sources of CP violation : 6l and the Majorana phases a, (3. 

To illustrate how these new contributions to Majorana phases occur, it is useful 
to recall again how we perform our calculations. Proceeding like in Section 6, using 



eqns. (117) and (119), we compute, from (38)(39), m2 and ms (with the conven- 



tion (pM|)) in terms of the given values for mi and tan6'i2. Then, m2 and 7713 get 
by construction new CP-violating contributions to the Majorana phases because. 



according to (117) (119), the inhomogeneous terms must be pure imaginary. 

Of course, since now we have new sources of CP violation in the Majorana phases. 



in the limit 5/, — we do not recover the simple limit (108) that we got for real 
values of 6*^3 and 6*^3 or, equivalently, for CP violation in the Majorana phases fixed 
exclusively from the phase 5l- 

We had CP violation in Majorana phases that were induced by their calculation 
for a given 5l in the case of vanishing C^g and C^^, (Subsection 5.2). But, from the 
imaginary inhomogeneous terms of the present Section, we have now new sources of 
CP violation in these phases. 

These new sources of CP violation in the Majorana phases, although small, are 



very efficient in producing a baryon asymmetry, as we realize from the results (118) 



and (120). 



It can be easily understood that the constraints (117) and (119) imply new 
contributions to the baryon asymmetry. These equations mean that the entries 



^23 = s-iid are purely imaginary. This in turn implies, from (43), new CP 



violation contributions to the mass matrix Mr, providing, after its diagonalization, 
new contributions to ei and Yb^- This important phenomenon certainly deserves 
further investigation for general complex inhomogeneous terms, keeping however a 
compact A''^ spectrum. 
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The important conclusion of the calculations of the present Section is that the 
results presented in Sections 6 and 7 remain much more general than the very 
particular model that, for the sake of simplicity, was exposed there. Our scheme, 
although fine-tuned because we look for a compact heavy neutrino spectrum, allows 
for a wide range of parameters giving good results. 

10 Comments on the compact heavy neutrino spec- 
trum and on the level crossing region 



We now go back to the case that we have studied in a quantitative detail, namely 



eqs. (|98j)(l99j) with the conditions (|101[)(|102[). 

Before and around the level crossing region we have a rather or very compact 
heavy neutrino spectrum. For simplicity, we have assumed that the lightest heavy 
neutrino N^-^^ decays out-of-equilibrium and gives the main contribution to the im- 
portant quantities relevant for baryogenesis : e, m and Yq. 

In such a fine-tuned situation, this can seem rather artificial. Actually, one 
should consider the contributions of all three heavy neutrinos, and therefore the 
contributions of all the CP-violation parameters ei, €2 and £3, and the corresponding 
wash-out factors. Notice that there are studies in the literature that consider all 
these contributions. See, for example, the paper by E. Bertuzzo et al. [26] and 
also, in a qualitative way, the work by Akhmedov et al. [7J. However, to take into 
account the contributions of all three heavy neutrinos can present some subtleties. 
We have not dared for the moment to roughly add Yg^, Yb^ and Yb^. We simply 
expect that the possible contributions of all three heavy neutrinos will not strongly 
affect the results that we have found from N^^ . An argument given below supports 
this hypothesis. 

Let us now comment on the crossing region. As pointed out above, there is a 
level crossing in both cases : (1) tan^s < and (2) tan 6*3 > 0. This happens around 
-C4 = 6-3^3 ~ 3 X 10*^ GeV-^ for solution (1) and = -Cg ~ 5 x 10^ GeV~^ for 
solution (2). At some point in this region the heavy neutrinos N^i and Nji^ become 
degenerate. 
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But we want to make explicit more precisely what we understand by level cross- 
ing. What we mean is that the properties of the N^-^ neutrino before the level 
crossing become (up to signs) those of the Npt^ neutrino after the level crossing, and 
vice-versa, exchanging their effect on the absolute magnitude of the final quantities 
e, m and Yg. This can be seen easily by writing the effective Dirac neutrino mass 



that enters in formulas (51), (54), before and after the crossing region. 

To illustrate what happens, we will take as an example solution (1). Similar 
features appear for solution (2). To be definite, we consider N^-^ (the hghtest neu- 
trino) and Nji,j (the next-to-lightest neutrino) before the crossing region, as we have 



computed in Section 7. Applying for Nr^ naively the formulas ( |54| ), (57) and (59) 
making just the exchange Mi o M2, let us give for solution (1) the quantities mi, 
ei, Ybi and 1712, €2, Yb^ at one point before the crossing region, for example for 
the value C^^ = 10^ GeV^^ and at one point after the crossing, for example for 



Cj"3 = 10^ GeV ^. Let us recall that the Dirac matrix (87) is completely fixed, but 



the redefined matrix (51) changes from point to point because it depends on the 



diagonalization of the matrix M/j by ( 49 ) . 

One finds, before the crossing, for C^^ = 10^ GeV~^, the values for the heavy 
neutrino masses : 

Ml = 5.531 X 10^ GeV M2 = 1.017 x 10^° GeV M3 = 2.017 x 10^° GeV 

(121) 



and the hermitian matrix that enters in (57) and (54) : 



0.259 
18.176 + 0.089i 
-0.125 + 25.597i 



18.176 - 0.089Z 
3352.06 
-0.00006 + 4720.59i 



and one gets therefore : 
nil = 0.047 eV 
m2 = 329.586 eV 



ei = -2.749 X 10"'^ 
€2 = -1.425 X 10"^° 



-0.125 - 25.597i 
-0.00006 - 4720.59i 
6647.84 



GeV 



(122) 
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Yb, = 2.383 X 10" 

Yb^ = 4.245 X 10"^° (123) 



Notice here one point. The numbers obtained in (123) are very interesting in 



relation with the calculations done in Section 7 for the region before the level cross- 
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ing. We have assumed there that Yb is dominated by the contribution of the hghtest 
neutrino Yb^. We see indeed that, at least within these naive estimate, this is true 
as far as the consideration of the next-to-hghtest neutrino is concerned. 



Remember that we have adopted the level ordering convention (53), that applies 



also after the crossing : Nji^ is the lightest neutrino and N^^ the next-to-lightest 
neutrino. 



One finds, after the crossing, for example for C; 
neutrino masses : 



L 
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10^ GeV-\ the heavy 



Ml = 2.011x10^ GeV Ms = 5.531 x 10^ Gel/ M3 = 1.020 x 10" GeV (124) 



and the hermitian matrix that enters in (57) and (54) : 
/ 



V 



193.27 

-5.937 + 0.020i 0.342 
-0.00029 - 1376.69i -0.143 + 42.293i 



-5.937 - 0.020i -0.00029 + 1376.69i 
-0.143 - 42.293i 
9806.55 



GeV^ 
(125) 



and one gets therefore : 



mi = 96.125 eV 



ei = 8.023 X 10 
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-9.981 X 10"^*^ 



ms = 0.062 eV 62 = 3.694 x lO"'' Yb, = -2.312 x 10"" (126) 

The shifts in order of magnitude among the elements of the matrices before and 



after the crossing, (122) or (125), explain the strong differences (in magnitude and 



even in sign) of the relevant quantities in these two regions, (123) or (126). We 
observe a strong discontinuity for the lightest neutrino properties (and for the next- 
to-lightest ones) that happens going through the crossing region. Up to the sign 
of e and therefore of Yb, we see that after the crossing the lightest neutrino has 
very strong wash-out and very small |ei| and therefore I^bJ, and that the opposite 
is true for the next-to-lightest heavy neutrino Nb,- It is easy to examine this for 
the parameter rhi (i = 1, 2), just by inspection of the matrix elements (mjm£))jj in 



(122) and (125). For it is a little more involved, but can be seen also by looking 



at the squares of the matrix elements of rh~^m£,- 

Let us now comment on the change of sign of e, and Yb^ before and after the 



level crossing region (123), (126). For solution (1), that we discuss here, the sign of 
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the inhomogeneous term (101) — C^3 = Cj^j > gives < and Yb^ > before the 
crossing, and > and Yg. < after the crossing. We have reahzed that if one 



changes the sign of ( 101 ), i.e. — = C^g < 0, then one has the opposite : ej > and 
Yb- < before the crossing, and < and Yb- > after the crossing. Adopting this 
latter sign for C^g, C^g, nothing essential changes for the heavy neutrino spectrum 



and for Am^ and Am^. The same considerations apply to solution (2) using (102) 
and changing its sign. 

Our conclusion is that, provided Nji-^ and are close enough in mass, one can 
have the right order of magnitude and sign for Yb before and after the crossing. 
This happens only at the price of changing the sign of our single free real parameter 

/^L /^L 

'^23 ~ '^33- 

An interesting conclusion is that, after the level crossing, the second-to-lightest 
heavy neutrino Npt^ dominates. The possibility of next-to-lightest neutrino domi- 
nance has been extensively studied recently by S. Antush et al. |27j . 



11 Open problems within the present approach 

There are a number of problems to face and study within the present approach. 
Let us make an incomplete list : 

(i) There is the possibility that more than one heavy right-handed neutrino de- 
cays out of equilibrium, contributing to the leptogenesis, a point that, in particular, 
has been suggested rather clearly in ref. [7]. If we guess a temperature T ~ 10^^ GeV 
below which all heavy neutrinos decay out-of-equilibrium, then not only the lightest 
Nr-^ decays out-of equilibrium after the level crossing, but also Nn^ is in the same 
situation. On the other hand, the heavy neutrino spectrum being rather compact, 
the natural thing to do would be to consider the contributions to Yb of all three 
heavy neutrinos Nji-^, N^^ and Nji^, i.e. to compute ei, €2 and eg and the relevant 
washout factors. 

For the moment we just expect that the consideration of the three neutrinos will 
not spoil the good features of the calculations of the present paper, that take only 
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into account the lightest neutrino A'^/j^ before the crossing. 

We have given an argument going in this sense in Section 10 where we have seen 
that, before the level crossing, the contribution of the next-to-lightest neutrino 
is negligible compared to the one of the lightest one N^i-^. 

(ii) One should take into account also the level crossing region and therefore the 
finite width of the right-handed neutrinos, as well as the delicate question of their 
interference. These problems have been treated in great detail by A. Pilaftsis and 
collaborators [25] that, to be complete, need to be adopted within our approach. 

(iii) The flavor effects, thoroughly studied by A. Abada et al. [28] are also a 
delicate question to study in this region of compact right-handed neutrino spectrum, 
and this should also be performed. 

(iv) It would be worth to study the more general case for CP violation outlined 
in Section 9, and make a detailed scan of the results in the case of general complex 
inhomogeneous terms - or equivalently general light neutrino Majorana phases -, 
with the constraint of having a compact heavy neutrino spectrum. 



(v) It could be that the homogeneous equations (60) correspond to some sym- 
metry, the inhomogeneous terms (that we have introduced to get a large enough 
CP violation ei) being a breaking of this symmetry, a possibility that would be 
interesting. 



12 Conclusions 



Our demand of a compact spectrum, and of an approximate quark-lepton 
symmetry implying a hierarchical spectrum for the Dirac neutrino masses with a 
similar structure between and the CKM mixing matrix, brings to a scenario 
where the lepton asymmetry comes out naturally, producing the required order 
of magnitude for the baryon asymmetry Yb ~ 0(10^^*^). We have assumed and 
justified that Yb is dominated by the contribution of the lightest neutrino Nn^. 

In this way, not only one can get a good magnitude for Yb, but as a natural 
consequence there are also a number of other strong points in this approach. 
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We get two possible solutions with a normal hierarchical light neutrino mass 
spectrum and an absolute scale, i.e. the lightest neutrino mass mi must be non 
vanishing. 

The light neutrino squared mass differences Am^ and Am^ are very stable and 
consistent with the data. 

There are three CP- violating phases in the whole approach, the phase 5^ of the 
unitary matrix, and the light neutrino Majorana phases a and (3. We take Sl 
to be of the order of the Kobayashi-Maskawa phase Skm- 

We have thoroughly studied in a quantitative way a particular case in which all 
CP violating effects are computed in terms of 5l, in particular ei and Yg^. It is 
interesting that one can get a baryon asymmetry of the right order of magnitude 
taking 6l — Skm- Of course, this result is not obtained in the Standard Model, but in 
a New Physics scheme under particular assumptions : Baryogenesis via Leptogenesis, 
5*0(10) Grand Unification, approximate quark-lcpton symmetry and compact heavy 
Nr spectrum. In the limit Sl ^ 0, one gets indeed ei ^ and Yb^ — )■ 0. 

The Ue mass, bounded by tritium /? decay, is of the order of few times 10~^ eV. 

The sum J2i i^ui satisfies the cosmological bounds, with a value rather close to 
the present upper limits. 

Let us emphasize that 5l induces also small CP violating corrections to the hght 
neutrino Majorana phases, that turn out to be naturally close to a = |, /? = | or 0. 
The effective neutrino mass, relevant for neutrinoless double beta decay, comes out 
to be rather small, of the order of 10~^ eV , because of strong cancellations due to 
the Majorana phases. 

In the region of quasi-degeneracy, the heaviest Nr has a mass of the order 1.5 x 
10^° GeV, roughly consistent with the expected scale ol B — L symmetry breaking, 
so that 50(10) breaks down to the Pati-Salam group 5C/(4) x SU{2) x SU{2) at 
the expected natural intermediate scale. 

We expose also an example in which the phase of the Dirac neutrino mass matrix 
5l and the Majorana phases a, /? are independent, providing an efficient generation 
of baryon asymmetry. 



46 



Acknowledgements 

This work has been supported in part by the EU Contract No. MRTN-CT-2006- 
035482, FLAVIAnet. We are indebted to A. Abada and F.-X. Josse-Michaux for 
discussions in the early stage of this work. 



Appendix 



In this Appendix we demonstrate the approximate formula (15) for Am^ 



Am2 = |m3p-cos2^^, Imsp-sin^^, \m^\ 



(A.l) 



The mass eigenstates read : 



\vi > = 


Cs 




> - 


Ss 

V2 




> + 


V2 




Ss 




> + 


Cs 




> - 


Cs 


> = 


1 
V2 




> + 


1 
V2 


\lJr 


> 





Wt > 



(A.2) 



and therefore the //-neutrino state is given, in terms of the mass eigenstates : 



> = - t^l^i > + ^ > + :h > 



V2 



V2 



(A.3) 



that evolves in time according to : 



fit,-. 



'"") \ / lit,) 



^i-ip-i^)t 



V2 ^ 



2p 



3 2\j- -, 



(A.4) 



and, from (A.2), we can write the scalar products : 



V2 



2 2 



2 2 

, — m2 



V2 



2 2 2 2' 

j(^o — -)i I «(— — 



2p 



V2 

2 2 

, m„ — m-, . 

1 + e ^ 2p ^ 



(A.5) 
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V2 



2 2 



s; e 



+ <e 



2p 



V2 
+ 1 



(A.6) 



V2 ^ 



2p 



2 



2p 



2 2 

c2 e^(^5^)* 



2i^i*+ 1 



Denoting 



q; 



(m2-m?)t 
2p 



(A.7) 



(A.8) 



One obtains, from (A.7) and (A.8) and the relation 



"31 -"32 = "21 



(A.9) 



the foUowing probabihties : 



< JyeWt^{t) > P = S^cl [1 - COs(a2l)] 



(A.IO) 



< u^Mt) > |2 ^ l[si + ci + l + 2s1c1 cos(q;2i) + 25^ cos(q;3i) + 2c1 cos(q;32)] (A.U) 

< u^Mt) > p ^ list + 4 + 1 + 2slcl cos(o;2i) - cos(q;3i) - 2c^ cos(a;32)] (A.U) 



and one obtains, as expected : 

I < UeWf^it) > P + I < l^nMt) > P + I < J^rMt) > P = 1 



— 2 



(m'2—m'f)t 
2^ 



Performing an expansion in powers of aji 

, , . -2_2 r^™2 ™2^+^2 

< ly^Mt) > |2 ~ 1 



2p 



', one finds 



4 



— 4 



2p 



2p 



+ - 



i r (m|-TOf)t 

1 [ 2p 

{m2—m\)t 
2p 



4 



(m|-m|)t 
2p 

(m|-m|)t 
2p 



(A13) 



(A.14) 

(A.15) 
(A.16) 
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with : 



The last terms in the r.h.s. of (A. 16) and (A. 17) reads : 



"2 r(m2— mj)t 
. 2^ 



2p 



and, for ml, ml << m^, the bracket in (A. 18) becomes 

[s^(m| — mf)^ + c^(m| — m2)^] ^ [7713 — (s^mf + 0^7722)]^ 
and therefore, formulas (A.14)-(A.16) can be approximated by 

I < Z/e|z/^(t) > |2 

I < u^Mt) > p 

I < lyrWi^it) > p 



{m'2—m'f)t 
2p 





" {m2—m'^)t' 


2 1 


"(m2-m2)t- 


4 


[ 2p J 


4 


L 2p J 



{m'^—m\)t 
2p 



+ i 



(m|-mg)f 
2p 



that satisfies (A. 17) and where 



777^ = slml+clm^ 



Therefore the formula (A.l) or (15) follows. 



(A.17) 



(A.18) 



(A.19) 



(A.20) 

(A.21) 
(A.22) 



(A.23) 
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